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ABSTRACT 

If a magnetic field is frozen into a plasma that undergoes spherical compression then 
the magnetic field B varies with the plasma density p according to B oc f?^^ . In the 
gravitational collapse of cosmological density perturbations, however, quasi- spherical 
evolution is very unlikely. In anisotropic collapses the magnetic field can be a much 
' steeper function of gas density than in the isotropic case. We investigate the distribu- 

, tion of amplifications in realistic gravitational collapses from Gaussian initial fluctu- 

ations using the Zel'dovich approximation. Representing our results using a relation 
of the form B oc p", we show that the median value of a can be much larger than 
' the a = 2/3 resulting from spherical collapse, even if there is no initial correlation 

lO . between magnetic field and principal collapse directions. These analytic arguments go 

' some way towards understanding the results of numerical simulations. 
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j>I ; 1 INTRODUCTION 
•I— I , 

. One of the most significant gaps in our understanding of cosmological structure formation is how large-scale magnetic fields 
?H ' arise and how they relate to the processes involved in galaxy and cluster formation. The observed galactic fields, with typical 
5^ 1 strengths of a few /xG, could be the result of the amplification of a tiny seed field, perhaps small as ~ 10~'^''G. On the other 
hand the adiabatic compression of a somewhat larger primordial field of ~ 10~^G could also achieve similar strengths in 
bound objects without violating cosmological constraints. For general reviews of various aspects of cosmological magnetic 
fields, see Kronberg (1994), Widrow (2002) and Vallee (2004). 

One specific situation where magnetic fields might prove important is in galaxy clusters. The most convincing line of 
evidence demonstrating the existence of magnetic fields in clusters emerges from studies of Faraday rotation measures (RM). 
Typical RM values of galaxy clusters are of the order of a few 100 rad/m'^, roughly consistent with field strengths of a few 
^G. These are comparable with galactic magnetic fields, but well below equipartition with the thermal cluster gas. Much 
larger RM values of a few 1000 rad/m^ have been detected in cooling-flow clusters, indicating possibly substantial magnetic 
pressure support of the intra-cluster gas there (e.g. Carilli & Taylor 2002) . Although the magnetic fields of galaxy clusters are 
less ordered than those of spiral galaxies, the presence of coherent structures is suggested by high resolution Faraday maps 
(e.g. Dreher et al. 1987; Taylor & Perley 1993; Taylor et al. 2001; Eilek & Owen 2002). These structures could result from 
shear-amplification of originally small-scale magnetic fields, as is suggested by MHD simulations of galaxy cluster formation 
(Dolag et al. 1999). It is clear, however, that the overall power-spectrum of magnetic field fiuctuations in clusters is rather 
broad ; see Dolag et al. (2002). The observed magnetic fields could play a significant role in supporting the intracluster medium 
through magnetic pressure (Loeb & Mao 1994) and there is some evidence that fields may be dynamically important in some 
clusters (Dolag et al. 2001; Eilek & Owen 2002). It is not known, however, how important these are or how they got there 
(Dolag et al. 1999; Goncavles & Friaca 1999; Dolag & Schindler 2000; Vlahos et al. 2005). 

In this paper we look at the behaviour of a primordial magnetic field during cosmological evolution. Our aim is to develop 
an understanding of how gravitational collapse can, even in the absence of turbulence or other dynamo action, lead to a 
significant amplification of magnetic fields during the the evolution of clusters or other large-scale cosmic structures. The 
essence of our argument is very simple. Assume there is a primordial magnetic field in a region undergoing compression 
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Figure 1. The deformation of a 2-D fluid clement of initial side lengths dqi, described by the deformation tensor TDij = The 
diagonal elements of the deformation tensor, & X'22 describe the scaling of the fluid element, while the off diagonal elements, 1512 
& X>2i describe the skew movement. 



and that the field is frozen into tlie plasma in tliat region. Since flux is conserved, the magnetic field strength scales with 
the cross-sectional area of the collapsing region while the gas density scales with its volume. The net result for a spherical 
compression is that B oc p^^"^ assuming no back- reaction of the field on the collapse (which we assume throughout this paper). 
One is tempted to infer that the "average" amplification due to collapse of a random set of density perturbations should be 
equal to this value. However, it has been known since the pioneering work of Zel'dovich (1970) and Doroshkevich (1970) that 
the generic collapse of cosmological perturbations is very anisotropic; see also Shandarin & Zelovich (1989). It is also true that 
the amplification of the B-field is a non-linear function of the overall collapse factor. Put these two facts together with the 
realisation that averaging does not commute with non-linear operations and it is clear that the angle- averaged amplification 
of B need not be close to the isotropic value. This argument was first presented by Bruni, Maartens & Tsagas (2003) who 
showed that there are possible collapse geometries that lead to a much larger boost in the B-field than one would naively 
expect. There is some evidence that this may be observed in galaxy clusters through a large value of the slope of the B — 1%^ 
relation: taking a model of the form B oc it seems values of a ~ 0.9 can be inferred (Dolag et al. 2001). The particular issue 
we are interested in tackling is what one can say about probabilities of these collapses and whether they lead to a significant 
effect after averaging over all possible collapse directions. 

Of course there are a number of complicated astrophysical processes that could boost the magnetic field in a cluster: 
turbulent dynamo action and the expulsion of galactic _B-fields are just two examples. These are beyond the scope of existing 
numerical methods, as well as analytic modelling, but it is nevertheless important to understand how much of the observed 
behaviour of cluster magnetic fields can be explained using simple arguments like those we present here. 

In the next Section we run through some basic definitions and simple examples of anisotropic collapse. In Section 3 we 
use an argument based on the Zel'dovich approximation to calculate the distribution of a for freely collapsing fiuid elements 
arising from Gaussian initial fiuctuations. 



2 BASICS 

Consider the effect of a general deformation of a fluid element. Define a deformation tensor, Vij, representing the deformation 
of an infinitesimal fiuid element of initial side lengths dqi, such that 

The deformation described by Dij is illustrated for the 2-D case in Figure Q 
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Figure 2. The effect of the deformation of a fluid element on the magnetic field lines in 2-D, assuming the field lines arc frozen into the 
fluid. The scaling of the fluid element affects the field strength only if the movement is perpendicular to the direction of the field lines. 
Skew movement affects both the field strength and direction. 



The initial co-moving volume of the fluid element, in three dimensions, is just dqidq2dq3, and its final comoving volume is 
\I)ij\dqidq2dq3, giving a ratio of final to initial volume of 

T T-fin ini 

^ = ^-12^-1- (2) 

The strength of a magnetic field is proportional to the density of its field lines. Hence, under the assumption that field 
lines are frozen into the fluid (as is usually assumed to be the case in astrophysical situations, due to the large scales involved), 
movement of the fluid perpendicular to the field lines affects the field strength, while movement parallel to the field lines leaves 
it unchanged. 



In order to see how the deformation tensor relates to the change in field strength it is easiest first to consider the two 
dimensional case, as illustrated in Figure |21 We first take the case where the initial field S'°' = only. The strength of the 
magnetic field is proportional to the area (or in this 2-D case, length) perpendicular to the field lines. Hence in Figure |51 

Srrfgi =5*""^/. (3) 

We can find Lf, the final length perpendicular to the field lines, by noting that 

Lf = asin(6l2 — 9i) = a (sin(S2) cos(Si) — cos(6'2) sin(Si)) , (4) 

so from Equation |3] the magnitude of the final field, B^", is given by 

ofln B'rdqi b BT 



a (sin(62) cos(6'i) — cos(6l2) sin(6'i)) dq2 



(5) 



9i = sm — cos I (6) 



V " 

„ . -1 (T)22dq2\ -1 f Vi2dq2\ 

92=sm (^^j=cos [^), (7) 

ICI = (©11022 -r'?2) • (8) 

Since the final field is no longer in the X2 direction only we also need to know its xi and X2 components, which are given by 

B^ = B«"cos(02) = Br^ (9) 



and 



and in two dimensions 
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and 

i3f = i3«"sm(e2) = Br^. (10) 

Following the same route for a field which is initially in the x\ direction only gives similar results, leading to total xi and X2 
components of the final field configuration given by 

Extending this argument to three dimensions gives 

5^ = ^E(sr^'»^)- (12) 



The total magnitude of a magnetic field is given by 



\B\ = \/(Si)' + (-B2)2 + (B3)2 (13) 
so in three dimensions the ratio of the magnetic field strengths before and after deformation is given by 



1 /E,,(Br)'^: 



(14) 

In the case of isotropic collapse, T>xx = T^yy — T>zz = 'Da, and all the ofl^ diagonal elements with i 7^ j are zero, giving 

finx 2/3 

(15) 



I _ '^ii I P 



(A,)' " Vp'"'. 

as expected, but the general case, which is not necessarily isotropic and can also include skew components, is more complicated 
than this. 

This is all very well, but we need to set up a model that applies to collapses of fluid elements within a spatially varying 
random fleld of density perturbations. Let us assume, for a toy model, that each possible direction of collapse is independent and 
equally likely. It is then straightforward to work out the average ampliflcation of the magnetic field which these deformations 
would produce. For reasons outlined in the introduction we focus on the parameter a, where 

Bocp". (16) 

We look first at the case of simple one-dimensional collapse, without skew components. If we assume that our initial magnetic 
field is in the x\ direction only, then for collapse in the x\ direction we have Dw — T>ii. The other diagonal elements 
2?22 = 2?33 = 1, with all off-diagonal elements 'Dij{i ^ j) ~ 0, which gives 



^2 



so that _B oc p" in this case. For collapse in the X2 or X3 direction only we have, for each case, 

IB**"! 1 



(18) 



In this case we see that B cc 1/V <x .We will from now on work with averages obtained over all possible configurations of the 
deformation by associating a value of a with each possible geometry and averaging these values. If these three eventualities 
discussed above are equally likely then the average value of a = 2/3, so that the average effect is B oc p^''^, just as in the 
isotropic case. 

We now look at the three possible skew directions, where the diagonal elements, 'Du = 1, and one of the off-diagonal 
elements T?ij = l^ij, with the other off-diagonal elements = 0. When is the only non-zero off-diagonal element we have 



_ 1 ViBiy + {Bl)HV,2r _ 1 

1+I?i2 vTbIF (1 + I?12)V2' ^ ' 

so that B oc p^^^. Similarly, when V13 is the only non-zero off diagonal element we also find B cx p^^^, but for the case where 
023 is non-zero, we find 

I D fin I 1 

f^ = 7-— , Bocp\ (20) 

(l+©23) ^ ' 

Again we find that the average over these three possibilities is B oc p^'''^, as for the case of isotropic collapse. 

However, in reality, each of these possible directions of collapse is not independent and equally likely. On the contrary, 
the elements of the deformation tensor are correlated, and as a consequence the amplification of the field is, on average, quite 
different from what would be expected from this naive model; see Bruni et al. (2003). In the next Section we investigate more 
realistic collapse geometries. 
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3 REALISTIC COLLAPSE MODELS 
3.1 The Zel'dovich Approximation 

The Zel'dovich approximation (Zel'dovich 1970; Doroshkevich 1970; Shandarin & Zel'dovich 1989; Bartelmann & Schneider 
1992; Salmi & Coles 1995) is an approach to structure formation which assumes that each particle has some small initial peculiar 
velocity which is determined by the initial density fluctuations. Particles then travel on ballistic trajectories determined by this 
initial veloicty, but unaffected by subsequent changes to the gravitational potential. Hence we write r{t), the final (Eulerian) 
position of a particle at time t, as 

r = a{t)q + a{t)b{t)v{q), (21) 

where a{t) is the scale factor of the universe, q is the initial (Lagrangian) position of the particle, b{t) is the growing mode, 
which scales the peculiar velocity, v, with time, and v depends only the initial position q. We can see that this equation looks 
like the background Hubble expansion, a{t)q, plus some perturbation which vanishes as t Q. 
Changing to co-moving co-ordinates, x — r/a{t), gives 

X ^ q + b{t)v{q). (22) 

If we assume that the initial velocity is irrotational (a reasonable assumption since it is generated by the gravitational effect 
of overdensities) we can write equation 12211 in terms of a potential, 0, giving 

x^q + b{t)\7<P{q), (23) 

where \/(j> = v, and V denotes differentiation with respect to the initial coordinates, q. This equation looks very much like 
ballistic motion, 

x = q + TV<P{q), (24) 

but with a peculiar time co-ordinate r — b(t). 

The Zel'dovich deformation tensor, "Dij, describes the deformation of a fluid element of initial side lengths dqi, and is 
given by 

T>ij = =5ij + T , (25) 

dqj Oqidqj 

where 5ij is the Kronecker delta and i,j = 1, 2, 3. The diagonal elements relate to the expansion and contraction of the fluid 
element along the principle axes, while the off-diagonal elements describe the skew distortion. The standard approach involves 
diagonalising the deformation tensor (which is symmetric). Let us denote the eigenvalues of ©ij by Afe, so that the mass 
density p(x, t) becomes 

p(x,f) =po(t)||^| =po(t)H?=i[l + 6(t)A,(q)]-\ (26) 

where po{t) is the mean cosmological background density. 

The Zel'dovich approximation describes the behaviour of evolving cosmological perturbations until the Jacobian in Equa- 
tion 12611 becomes zero, at which point the matter density becomes infinite and a caustic forms, where the mapping between 
Eulerian and Lagrangian coordinates is no longer unique. Although simple, the Zel'dovich approximation gives remarkably 
accurate results up to shell crossing (Coles et al. 1993). This, along with its natural formulation in terms of the deformation 
tensor describing the collapse of each fluid element, makes it an ideal approach for studying the effect of gravitational collapse 
on magnetic fleld strength. 



3.2 Probabilities 

The joint probability density function (JPDF) for the elements of the Zel'dovich deformation tensor is given by Doroshkevich 
(1970), where a very elegant discussion is also presented about the probability of various collapse geometries. The most likely 
collapse is initially along one direction, resulting in caustics that are generically two-dimensional ("pancakes"). To proceed 
we assume the initial density perturbations are Gaussian, as is expected from inflationary models (Bardeen et al. 1986). The 
initial velocity potential is then also Gaussian, and so are its derivatives. Therefore the JPDF of the second derivatives of the 
velocity potential, on which the elements of the deformation tensor depends, is a multivariate Gaussian distribution. If we 
write all the possible second derivatives as a vector, ipp, p = 1...6, with the diagonal elements, which have i = j, as (pi, ifi2, <P3 
and the off-diagonal elements with i 7^ j as ^4, ^5, ^6, then their joint probability is given by 

|C-/|5 exp -i J2p,q ^pCpqV>g\ 

P(V.1...<^6) = ^ (27) 

(27r) 
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and we need only find the value of the covariance matrix, Cpg, to know the probability distribution in full. Note that all the 
variates involved have zero mean. The components of Cpq are given by 

^ {^pipq) = (^-^^-^> = ^ + S.,S,, + SuS,,) , (28) 

= f P{k)k'dk, (29) 



where 

2 

rr = 

27r2 _ 

if P{k) is the power spectrum. So we now know the multivariate Gaussian distribution which gives the JPDF of the elements 
of the deformation tensor, provided we know P{k), the power spectrum of the initial density perturbations. For simplicity we 
take the case where the magnetic field is entirely independent of the deformation tensor (i.e. we are not appealing to the field 
to create the initial deformation in any way, so there is no correlation between the field and the deformation) , and ignore any 
back reaction of the field on the fiuid. This gives us the freedom to choose to have the initial magnetic field along one axis, 
eg B'"' — BJ"' only. This simplifies the equation for the field strength amplification (Equation considerably, giving 

Since the marginal distribution P(x) — J P(x\y)P(y)dy, the probability distribution of the magnetic field amplification, given 
the deformation tensor, is given by 

Pff l-'""!^ /°°A ( 1 ^7^2 ^T)2 ^T^M^/^ |C-/|^Cxp[-iyJC-,V,] 



where 



l^-l = ^fS). (32) 



and Soi^:) is the Dirac delta-function. Similarly, the the density compression factor is 



„fln 



/p = P4 (33) 



and its probability distribution is 



|C-/|^exp[-ly>JC-,\ 



The two distributions (I31II and (I34II are both skewed showing that there are amplifications and compressions much larger 
than average. This skewness also increases with time, as the perturbations become more non-linear. This suggests that there 
is a significant probability of large amplifications of B relative to p, as suggested by Bruni et al. (2003). However these two 
distributions are not sufficient to calculate the distribution of our chosen parameter a, which we define by 

" - (35) 

which requires the full joint distribution P{fB,fp)- Unfortunately the integrals involved in the construction of this joint 
distribution are resistant to analytical methods. In order to investigate the statistical properties of a (where B cc p"), we 
therefore resorted to numerical techniques. 



3.3 Calculations 

Although the multivariate distributions described above are difficult to handle analytically, it is fortunately quite easy to obtain 
the distribution of a using a Monte Carlo technique. First we need to generate deformation tensors from the distribution 12711 . 
One can generate correlated Gaussian variates with a given covariance matrix by simply reversing the process by which one 
can diagonalise the covariance matrix to produce independent variates from correlated ones. In other words, we start with a 
set of six independent and identically distributed Gaussian variables, and then construct a linear combination such that they 
produce the correct covariance matrix. This is easily done using a Cholesky decomposition. We use this technique to generate 
a large library of deformation tensors that constitute a fair sample from the required distribution 1271 . We then used the 
machinery described above to calculate the compression factors /s and fp from this sample. This in turn generates a value 
of a for each fluid element from which the overall distribution can be formed. For simplicity we have assumed, as before, that 
there is no correlation between the magnetic field and the deformation, and have ignored any back reaction of the magnetic 
field on the fluid. 

There are a number of subtleties that must be confronted before dealing with this distribution. The flrst is that not 
all deformation tensors compatible with our parent distribution lead to collapsing fluid elements. Since we wish to look 
at collapsed regions only, we consequently restrict our library of deformations describing a net compression in comoving 
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Figure 3. Frequency distribution of fs when <5 = 1. The long tail of the distribution (which continues to the right of what is shown 
above) is due to chance alignments of the direction of collapse with the magnetic field lines which have a low but non-zero probability. 



coordinates. Secondly, since our variates are Gaussian the eigenvalues Xi are also Gaussian, Equation (I26II always assigns a 
non-zero value to regions where shell-crossing has occurred. This generates a tail in the distribution that is unphysical. For 
this reason we illustrate our results by focussing on the field strength amplification which would be generated by a freely 
collapsing fluid element whose overdensity, 5, has reached a value 



S = 



(P) 



(P) 



= 1. 



(36) 



Each fluid element in our library was allowed to collapse to this state, at which point the magnetic held compression factor, 
/s, which depends strongly upon the exact conflguration of the collapse, was determined using methods given in Section 2. 
This was then compared to the change in density to flnd a value of a for each element. These values were binned to give 
flnal distributions of /b and a which can be seen in Figures |3 and |1| respectively. Notice that the distributions of both B 
and a have a very long tail of values. This effect is understandable in terms of a simple model: if a fluid element collapses in 
directions orthogonal to the magnetic field, the field obviously undergoes a large boost. However, it is possible for the fluid 
element to be expanding along the third direction (along the field). This means that B can increase significantly while p 
hardly changes at all. This can lead to unbounded values of a. These are unlikely, but not impossible so they have a finite 
probability. The mean value of a is consequently formally undefined. 

Intriguingly, the median value of a once all the fluid elements had collapsed until 5 = 1, is found at a value of a = 1.19, 
considerably higher than the expected isotropic value of a = 0.67. One has to be careful interpreting this result as it does not 
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Figure 4. Frequency distribution of a when 5 = 1. Tlie mean value of a is undefined due to tiie long tail in the distribution, which is 
caused by chance alignments of the direction of collapse with the magnetic field lines. The median value of a is 1.19. 



correspond to a calculation of a at the same epoch. Some fluid elements will reach 5 = 1 at different times to others. Since 
the background B oc and the background p (x the initial ratio effectively depends on the time taken to collapse 

1/3 

through a{t) or, equivalently as pg (t). To compare these results at the same epoch one should therefore subtract 1/3 from 
the result above giving a ~ 0.86. 

While most deformations cause a fleld ampliflcation at or near the median value (approximately 87% are within 50% of 
the median), there are some outliers caused by unusual alignments of the magnetic field compared to the direction of greatest 
collapse as discussed above. High values of a are caused by configurations where the collapse of the deformation tensor is 
such that the magnetic field is amplified greatly for only a small change in the density. This occurs when the magnetic field is 
perpendicular to the axis of fastest collapse, and parallel to an axis along which the collapse is slow, or expansion is occurring. 
On the other hand, if this happens the magnetic field may reach such a large level that it begins to exert a back-reaction on 
the collapse. Our treatment can not handle this eventuality accurately. 

Similarly, extremely low values of a, some of which are negative, occur when the field is perpendicular to an axis which 
is collapsing slowly, or even expanding, and parallel to a quickly collapsing axis, causing a large increase in density for very 
small increase, or even a decrease, in magnetic field strength. Such alignments are unusual, but do occasionally occur, and 
cause the long tails in the distribution of a seen in Figure 2] 
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4 DISCUSSION 

We have investigated the probable magnetic field amplification by anisotropic gravitational collapse under a number of 
restrictive assumptions. These are: 

(i) the initial fiuid perturbations are Gaussian; 

(ii) the initial magnetic field is uncorrelated with the perturbations; 

(iii) collapse is described by the (quasi-linear) Zel'dovich approximation; 

(iv) shell-crossing does not occur; 

(v) there is no back-reaction of the magnetic field on the fluid motion. 

We find the median value of a to be significantly higher than the isotropic value of a = 2/3 — 0.67, but the precise value 
depends on exactly how one performs averages. We have taken two definitions (which differ by 1/3 in q) but other ways of 
averaging are possible. Indeed, for comparison with a single cosmological object one would probably have to average over 
correlated fiuid elements, and average by mass rather than by direction. Obviously if one averages over coherent structures 
and/or if there are correlations between magnetic field and the fiuid perturbations, such as would be the case if the field were 
generated by fiuid activity, then the results would be very different. 

In the numerical simulations presented by Dolag et al. (2001) it seems values of a ~ 0.9 can be reached. We cannot 
argue that our simplified treatment can be compared directly with these more accurate computations, but it is reassuring that 
taking proper account of collapse geometries leads to plausible values of a that are not far from those reached in fully MHD 
approaches. This suggests at any rate that it may be possible to get some insight into the behaviour of cluster magnetic fields 
using simple arguments like those we have presented here, without needing to appeal to more complex mechanisms such as 
sheer flows and dynamos. 
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